Quantum transport through a DNA wire in a dissipative environment 
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Electronic transport through DNA wires in the presence of a strong dissipative environment is in- 
vestigated. We show that new bath-induced electronic states are formed within the bandgap. These 
states show up in the linear conductance spectrum as a temperature dependent background and lead 
to a crossover from tunneling to thermal activated behavior with increasing temperature. Depend- 
ing on the strength of the electron-bath coupling, the conductance at the Fermi level can show a 
weak exponential or even an algebraic length dependence. Our results suggest a new environmental- 
induced transport mechanism. This might be relevant for the understanding of molecular conduction 
experiments in liquid solution, like those recently performed on poly(GC) oligomers in a water buffer 
(B. Xu et al, Nano Lett. 4, 1105 (2004)). 
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In the emerging field of molecular electronics, DNA 
oligomers have drawn in the last decade the attention 
of both experimentalists and theoreticians [1], This has 
been mainly motivated by DNA exciting potential appli- 
cations which include its use as a template in molecu- 
lar devices or by exploiting its self-assembling and self- 
recognition properties [2]. Alternatively, DNA strands 
might act as molecular wires either in periodic confor- 
mations as in poly(GC), or by doping with metal cations 
as is the case of M-DNA [3] . As a consequence, the identi- 
fication of the relevant charge transport channels in DNA 
systems becomes a crucial issue. Transport experiments 
in DNA derivatives are however quite controversial [4, 5]. 
DNA has been characterized as insulating [6], semicon- 
ducting [7] or metallic [8, 9]. It becomes then appar- 
ent that sample preparation and experimental conditions 
are more critical than in transport experiments on other 
nanoscale systems. Meanwhile, a variety of factors that 
appreciably control charge propagation along the double 
helix have been theoretically identified: static [10] and 
dynamical [11] disorder related to random base pair se- 
quences and structural fluctuations, respectively, as well 
as environmental effects associated with correlated fluc- 
tuations of counterions [12] or with the formation of lo- 
calised states within the bandgap [4, 13]. 

Recently, Xu et al. [9] have carried out transport ex- 
periments on poly(GC) oligomers in aqueous solution. 
These experiments are remarkable for different reasons: 
(i) it was shown that transport characteristics of sin- 
gle molecules were probed, (ii) the molecules displayed 
ohmic-like behavior in the low-bias I-V characteristics 
and (iii) the linear conductance showed an algebraic de- 
pendence g ~ iV _1 on the number N of base pairs. This 
latter result suggests the dominance of incoherent charge 
transport processes. Complex band structure calcula- 
tions [14] for dry poly(GC) oligomers predict, on the con- 
trary, a rather strong exponential dependence of the con- 
ductance on the wire length, a typical result for coherent 
tunneling through band gaps. Hence, Xu et al. exper- 



iments suggest the potential role of the environment in 
modifying the DC conductivity of DNA. 
In the light of these results, we explore in this Letter 
the possibility that a strong perturbation of the elec- 
tronic system by a dissipative environment may lead to a 
modification of the low-energy electronic structure of the 
molecular wire. As a result, the linear transport proper- 
ties may be qualitatively different when comparing with 
the "dry" wire. 

Our description assumes that only the frontier orbitals 
of the poly(GC) stack are relevant for charge transport, 
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FIG. 1: Schematic representation of the DNA molecular wire 
in contact with a phonon bath (upper panel) and of the corre- 
sponding density of states (lower panel). In the absence of the 
phonon bath, valence (VM) and conduction (CM) manifolds 
are separated by a gap. Upon coupling to the environment, 
a new set of states emerge within the gap. Though strongly 
damped by the dissipative coupling, they yield a finite density 
of states and eventually lead to a conductance enhancement 
with increasing temperature. 
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a reasonable approximation at low bias. Frontier or- 
bitals are the highest-occupied (HOMO) and the lowest- 
unoccupied (LUMO) molecular states. They both have 
7r character and are derived from linear combinations 
of the p z orbitals of individual base pairs. The 7r or- 
bital stack can be thus represented in a localised or- 
bital picture. As shown by first principle calculations 
[15], the HOMO charge density is mainly concentrated on 
the guanine bases and the LUMO charge density on the 
citosyne bases. Hence, within a minimal model, we will 
focus e.g. only on the 7r-stack along the guanine strand 
(HOMO) for hole transport and consider the opening of 
the HOMO-LUMO gap as a perturbation of this 7r-stack, 
induced by the complementary strand and, eventually by 
the backbone subsystem. [] 6] The environment is conven- 
tionally described by a harmonic phonon bath. 
We address the temperature and length dependence of 
the conductance in the strong coupling limit to the bath 
degrees of freedom. Our results can be summarized as 
follows. First, bath-induced states appear in the gap re- 
gion, see Fig. 1 for illustration. They are however washed 
out due to the strong effect of the environment, so that 
they do not manifest as well-defined resonances in the 
transmission spectrum. Nevertheless, they induce a tem- 
perature dependent density of states within the gap. This 
leads to a transition from a tunneling regime at low tem- 
peratures, with a zero current gap, to an activated regime 
at higher temperatures, with a finite current at low volt- 
ages. Second, a weak exponential or even algebraic length 
dependence together with an Arrhenius-likc behavior of 
the transmission at the Fermi energy are found, reflecting 
the strong contribution of incoherent processes. 
We describe the system consisting of a poly(GC) wire 
containing N base pairs, contacted to left and right elec- 
trodes (7iieads); and in interaction with a phonon bath 
(Hb) by the following Hamilton operator, sec Fig. 1: 

T~l = 7~Lc + %c + Hloads + Wb, (1) 

■He = e b J2 % - t„ J2 ( 6 5 & i+i + Rc -) > 

3 3 

Hc-c = <X'V-< -*iE t'V' + "•'••) > 

3 3 

Pleads = e ko- rf ka rf ko- + ^ f^M & 1 + H.e.J , 

keL,R,cr k6L,cr 

+ ( Vk > N d ka b N + H.c.) , 

keR,cr 

In the previous expression, TCc describes the HOMO ir- 
stack as a one-dimensional tight-binding chain with on- 
site energies e b and intra-strand hopping in. The per- 
turbation arising from the complementary strand is de- 
scribed via the side-chain Hamiltonian TLc-c- The inter- 
chain hopping tj_ is assumed to be a small parameter, 
according to recent numerical estimates. [] 7] As a result, 



the probability of inter-chain charge hopping will be very 
small and we may thus neglect in a first approximation 
charge propagation on the side chain, i.e. charge trans- 
port occurs only on pathways along the central chain. 

In Eq. (1), the onsite energies e& and e will be set equal 
to zero for a homogeneous wire. When simulating the 
introduction of A— T bases on a finite segment of the 
molecular wire, we set e b ^ 0. The set of bath fre- 
quencies fl a and corresponding coupling constants A Q , 
q : = 1, . . . , M, does not need to be further specified. By 
performing the thermodynamic limit (M — ► oo) later 
on, the bath can be described by a spectral density 
A") = E a \ 2 J(Lu-n a ) = J (u/LJ c ye-"/"°G(u), where 
lu c is a cut-off frequency and Q(u>) is the Heaviside func- 
tion. [18] In what follows, we consider only the case s = 1 
which corresponds to an ohmic bath. The form of the 
spectral density will of course depend on the specific type 
of environment. In the case of an aqueous environment, 
a Dcbye model for dielectric relaxation might seem more 
appropriate. [19] However, its low-frequency behavior is 
similar to that of the ohmic bath; hence, we can safely 
approximate it by the above spectral density. We do not 
expect that the exponential cut-off at high frequencies 
will have a very dramatic effect on our results. 

For X a = Va, the model shows a temperature indepen- 
dent gap in the electronic spectrum, the gap being basi- 
cally proportional to t±. Valence and conduction mani- 
folds, involving N states each, are symmetric w.r.t. the 
Fermi level which is the zero of energy (particle-hole sym- 
metry) . For nonzero coupling to the bath a more involved 
behavior may be expected depending on the electron- 
phonon coupling strength. The interaction with the bath 
degrees of freedom can be eliminated by performing a 
unitary transformation [18]. As a result the transversal 
coupling Hamiltonian Hc-c will be renormalizcd by the 
bath operators [20]. 

Since we are mainly interested in the temperature and 
length dependence of the linear conductance, we ad- 
dress only the low-bias regime, thus neglecting nonequi- 
librium effects as well as inelastic tunneling which be- 
come more relevant at large applied voltages. As a re- 
sult, the current can still be written as [21] I(V) = 
(2e/h) j dE (f(E - eV/2) - f(E + eV/2)) t(E). Note 
that the function t(E) = 4Tr {lmS L G ImS R } con- 
tains contributions from the electron-bath interaction via 
the wire Green function G. It hence describes processes 
where the incoming and outgoing electron energies are 
equal, though virtual phonon emission and absorption 
is allowed and included to all orders in G. We denote 
in what follows t(E) as a transmission, though it is not 
directly related to the usual Landaucr transmission func- 
tion due to the presence of the dissipative bath in the wire 
propagator. 

The Green function G can be calculated using equation 
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FIG. 2: Upper panel: The function t(E) for two different 
temperatures; the inset is a log-plot around E = showing 
the strong temperature dependence of the pseudo-gap. Lower 
panel: I-V characteristics. Parameters: N = 20, Jo/W = 
20, *, | =0.6eV,tx/<|| =0.2,r L/R /f N =0.16. 

of motion techniques. One finds to lowest order in t±: [20] 

G-^E) = El-H c -'E{E)-tlP(E), (2) 

P tj {E) = % / die 1 '^""^G^e "*<*>, 
Jo 

where £(£) = £ L (£) + £ R (£), e"*« = (#(t)#t(0)) B 
is a dynamical bath correlation function and X = 
cxp \Y^ a (\ a /£t a )(B a — £?£)] . The electrode selfcnergies 
Sl/r are calculated in the wide-band limit, Hi, i £j(E) = 
-iTx^Suhj and T lR j j (E) = -iT R S Ni d N j, i.e. ignoring 
their energy dependence. Note that the function P{E) 
containing the free electron Green function G° (t) of the 
side chain and the bath correlator exp (— $(t)), act as an 
additional selfcnergy for the central chain propagator. 
In the weak-coupling regime to the bath, no relevant 
physical effects were found [20]. In the following, we 
discuss the strong-coupling limit Jq/lo c > 1, where an 
appreciable modification of the electronic spectrum oc- 
curs. Recent estimates [22] of the latter parameter using 
the classical Onsager model for molecule-solvent interac- 
tions suggest that this regime can be realized in a water 
environment. 

In Fig. 2 the transmission t(E) and the corresponding 
current are shown. A strong temperature dependent gap 
in the electronic transmission spectrum t(E) is found, its 
magnitude increasing with temperature. The low-voltage 
I-V characteristics evolve from a "semiconducting" be- 
havior at low temperatures (zero current around V = 0) 
to a "metallic" behavior (nonzero slope near V = 0) 
with increasing temperature. The reason is that in the 
strong dissipativc regime a pseudo-gap rather than a gap 
in the electronic spectrum is induced by the bath dy- 
namics. An analysis of the real and imaginary parts of 
P(E), Eq. (2), at low energies helps to understand this. 
One can show that (i) ReP(E) ~ E for E ~ and (ii) 



ImP(E) is peaked at E = 0. For comparison, in the 
absence of the bath ReP(E) would display a 1/E be- 
havior around E = [16]. It follows from (i) that addi- 
tional low energy poles of the wire Green function G(E) 
might emerge symmetrically placed around the Fermi en- 
ergy, building a third electronic manifold. These states 
can show up as resonances in the transmission spectrum 
inside the gap, see Fig. 1, and form a polaronic band. 
We stress that they are neither present for X a = nor 
in the weak-coupling limit. We would then have three 
electronic manifolds. It turns out, however, that the 
non- vanishing Im P(E) (the "frictional" part) has a dra- 
matic influence on these states. Since they are located 
in an energy region where Im P(E) is appreciably dif- 
ferent from zero, no well-defined resonances manifest in 
the low-energy sector of the transmission. Nevertheless, 
these bath-induced states do contribute with a temper- 
ature dependent incoherent background and eventually 
lead to an increase in the density of states near E = 
when the temperature grows. Hence, the current may be 
enhanced at low voltages with increasing temperature. 
We thus interpret the central manifold as an incoherent 
polaronic band that supports activated hopping at high 
temperatures. This is farther reflected in the tempera- 
ture dependence of the transmission at the Fermi energy, 
which displays an Arrhenius-like law, see Fig. 3. We note 
in passing that a qualitative similar effect is found in the 
so called spin-polaron problem. [23] 
A controversial issue in transport through DNA-based 
systems is the actual length dependence of the electron 
transfer rates or correspondingly, of the linear conduc- 
tance [24, 25, 26]. Different functional dependences have 
been found in charge transfer experiments ranging from 
strong exponential behavior related to superexchange 
mediated electron transfer [25] to algebraic dependences 
typical of thermal activated hopping [24, 26]. As far as 
transport experiments are concerned, Xu et al. [9] re- 
ported an algebraic length dependence of the conduc- 
tance for poly(GC) oligomers in solution. Theory has 
shown that a transition between different regimes may 
happen as a function of the wire length [27]. We have 
investigated the length dependence of t(Ep) and found 
different scaling laws dependent on the strength of the 
electron-bath coupling. For Jq/lu c > 1, an exponen- 
tial law for energies close to E-p was found, t(E-p) ~ 
cxp(— 7L), see Fig. 3, middel panel. Here, L = Nclq, 
ciq ~ 3.4 A being the average inter-base separation. At 
the first sight, this might be not surprising since a gap 
in the spectrum does exist. Indeed, in the absence of the 
bath, i.e. with an intrinsic electronic gap, we get decay 
lengths 7 co h of the order of 2 A . However, as soon as 
the interaction with the bath is included, we find values of 
7 much smaller than expected for virtual tunneling, rang- 
ing from 0.15 A to 0.4 A . Additionally, 7 is strongly 
dependent on the strength of the electron-bath coupling 
as well as on temperature; 7 is reduced when Jq/lu c or 
I^bT increase, see Fig. 3, since in both cases the density 
of states within the pseudo-gap increases. These results 
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FIG. 3: Upper panel: Arrhenius plot for t(Ep). Parameters: 
TV = 20, f|| = 0.6eV, t±/t\\ = 0.2, F L/R /t u = 0.16. Inset: 
temperature dependence of the exponential decay length 7 
(see text). Middle and lower panels: Length dependence of 
t(Ep) at different temperatures for two different strengths of 
the electron-bath coupling Jo/uj c . The electronic coupling pa- 
rameters are the same as in the upper panel. The inset in the 
lower panel shows the effect of introducing tunnel barriers to 
simulate the insertion of Vat A-T base pairs in an otherwise 
homogeneous chain of length N (T=300 K). The algebraic 
lenght dependence goes over onto an exponential one. For 
this, the onsite energies along a finite segment of the wire 
were shifted by £b = — 1.5 eV, see Eq. (1). The number of 
unperturbed sites Ngc in a sequence like Agc — Vat — Vgc 
was kept constant (Ngc — 4) while varying Vvr = 1 ■ ■ ■ 8. 

clearly indicate that the bath does strongly determine 
the effective decay length in a way which we can quan- 
tify by extracting a 7 cn v term such that 7 = 7 co h — 7cnv- 
The first contribution 7 co h is purely determined by the 
intrinsic electronic structure of the wire and can be ob- 
tained, e.g. by means of complex band structure ap- 
proaches [14, 28]. A 7 co h of the order of 1.5 A has 
been recently calculated for poly(GC) [14], which com- 
pares well with our estimated 7 co h- The dependences on 



Jo and /ce? 1 are hence contained in the bath-induced con- 
tribution 7 cnv . Remarkably, in the regime Jq/lu c 1 the 
weak exponential length dependence goes over onto an 
algebraic dependence, see Fig. 3, lower panel. The in- 
troduction of a tunnel barrier induced by the insertion 
of A— T base pairs in the poly(GC) oligomer [9], can be 
simulated in the simplest way by a shift of the onsite en- 
ergies along a finite segment of the wire. In this case, 
the exponential dependence is recovered (see the inset of 
Fig. 3, lower panel). 

In conclusion, we have investigated the influence of a dis- 
sipative environment on charge transport along a molec- 
ular wire in a model that mimics basic features of the 
electronic structure of poly(GC) oligomers. We found 
a strong modification of the low-energy electronic struc- 
ture of the wire in the strong dissipative regime. Indeed, 
a pseudo-gap is formed which induces a temperature de- 
pendent background around the Fermi energy. The re- 
sulting non- vanishing low-bias current at room tempera- 
ture as well as the algebraic dependence of the conduc- 
tance on the wire length found in our model suggests 
that the striking results of Xu et al. transport experi- 
ments [9] may be related to the presence of an aqueous 
environment. The interplay with other dynamical de- 
grees of freedom like internal molecular vibrations has, 
of course, to be further clarified. We note that the in- 
clusion of randomness in the base pair distribution (as is 
the case of e.g. A-DNA) docs not qualitatively change the 
above picture [20]. Disorder mainly washes out the side 
bands in the transmission without essentially changing 
the behavior around the Fermi level. 

Finally, we remark that a close estimation of the 
physically relevant model parameters, especially of the 
electron-bath interaction strengths, requires a detailed 
analysis of first principle calculations of DNA oligomers 
in solution, eventually including vibrational degrees of 
freedom. This goes, however, not only beyond the scope 
of this investigation, but also beyond the actual capabil- 
ities of most ab initio approaches. 
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